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Abstract 

We observe the anisotropy of the power spectral tensor of magnetic field fluctuations in the 
fast solar wind for the first time. In heliocentric RTN coordinates the power in each element of the 
tensor has a unique dependence on the angle between the magnetic field and velocity of the solar 
wind (8b) and the angle of the vector in the plane perpendicular to the velocity (4>b)- We derive 
the geometrical effect of the high speed flow of the solar wind past the spacecraft on the power 
spectrum in the frame of the plasma P(k) to arrive at the observed power spectrum P(f, 0bi4>b) 
based on a scalar field description of turbulence theory. This allows us to predict the variation 
in the <\>b direction and compare it to the data. We then transform the observations from RTN 
coordinates to magnetic-field-aligned coordinates. The observed reduced power spectral tensor 
matches the theoretical predictions we derive in both RTN and field-aligned coordinates which 
means that the local magnetic field we calculate with wavelet envelope functions is an accurate 
representation of the physical axis of symmetry for the turbulence and implies that on average 
the turbulence is axi-symmetric. We also show that we can separate the dominant toroidal 
component of the turbulence from the smaller but significant poloidal component and that these 
have different power anisotropy. We also conclude that the magnetic helicity is anisotropic and 
mostly two-dimensional, arising from wavevectors largely confined to the plane perpendicular to 
B. 



1 Introduction 



Fast solar wind from the poles of the Sun is an excellent example of MHD turbulence, wi th the 
fluctuations being approximately incompressible ( Goldstein et al. . 1995 ; Horbury et al. . 20051 ) . The 
Ulysses spa cecraft provides a unique data set with extended periods i n this continu o us fa st polar 
solar wind ( Ebert et al. . 20091 ) and high cadence magnetic field data ( Balogh et al. . 1992 ). Such 
observations allow us to investigate how turbulence makes the nominally collisionless solar wind 
behave like a gas with shocks and structures, and why superthermal particles and cosmic rays 
appear to be diffusively coupled to the solar wind, allowing exchange of energy. The details of that 
coupling are not yet understood completely, and the poorly understood anisotropy of the turbulence 
is a part of the problem. Recently measurements have been made using th e Ulysses data c learly 



showing the importance of the magnetic field dire ction in the turbulence (jHorbury et all 120081 : 



Podestal . 120091 : ILuo fc Wul . l20ld : IWicks et all I2OI0I ) , with different power amplitudes and spectra] 
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indices in different directions relative to the local mean magnetic field. Attempts have been made 



to choose between theories of anisotro pic turbulent cascades (e.g. iGoldreich fe Sridharj . Il995l . 119971 ; 



BoldyrevL 120061 : iLithwick et al.l . 120071 ) by observing the scaling of the power at different angles to 



the mean field. 

All of these studies have concentrated on the trace of the magnetic power spectral tensor (i.e. the 
total power in magnetic fluctuations) rather than the whole tensor, which is needed to fully describe 
turbulence. The second order correlation tensor and the associated power sp e ctral t ensor are central 
parts of generalized turbulence theories (jRobertsonl . Il940i : iBatchelorl . Il946l . Il970l : IChandrasekharl . 
195Ch . Incompressible MHD turbulence is different in many respects to incompressible hydro d y- 



namic turbulen ce, primarily since it has two solenoidal fields, V and B ( Chandrasekhar , 1951al b: 



Biskampl . [2003h . Theoret cal treatments show that solenoidal fields in MHD plasmas (e.g. V--B = 0) 



require correlation and power spectral tenso rs which are complete ly described by four standard ten 



sor forms multiplying four scalar functions (jOughton et al.l . Il997l ) 



Frequency power spectra from single-spacecraft observations P(f) ar e equivalent to the 'reduced ' 
form of the full thr ee-dimensional wave vector power spectrum P(k) (jFredricks fc Coronitil . Il976l : 
Forman et all l201ll ). In this context, k represents the wavelength and orientation of the 3D spatial 



struc ture of the turbulence, which is advected past the spacecraft at supersonic speeds ( Tavloi . 
1938). The resulting integral is a type of tomographic projection called a Radon transform (jRadonl . 



19171 : iDebnatha fc Bhattal . 120071 ). It is impossible to separate wave vectors that h ave the same 
projection on the direction of flow ( Fredricks Sz Coroniti , 19761 : Forman et al. . 2011 ). This causes 
a permanent ambiguity in the observed power spectrum and means that in- situ observations by 
single spacecraft can only fully resolve the 3D spectrum or the related correlation tensor if they are 
isotr opic. Multiple-spacecraft missions, such as C l uster, have been used to overcome this problem 
(e.g. 



Osman fc Horburv , 2007 : Narita et al. . 2010l : Sahraoui et al. . 2010l ) but do not spend long 



m 



the solar wind and so are difficult to use for turbulence studies which require ensemble averages 
over large data sets, and the resolution of wave vectors is relatively coarse. 

In Section 2 we use the Ulysses data to make the first measurements of all nine elements of the 
reduced magnetic power spectral tensor Pij(f, b) as a function of the direction b of the local mean 
mag netic field in the solar wind. These are measured in the Sun-spacecraft aligned RTN coordinates 
(see iBurlagal . Il984l : iFranz fc Harperl . [2002J) and show a remarkable amount of variation with b. In 
Section 3 we derive the dependen ce of the reduced Pij(f,b) on Pij{k) using the general tensor 
formalism of lOughton et al.l ( 19971 ) and show how the resulting 4 scalar functions appear in the 
measured power spectral tensor in RTN coordinates. This may seem inelegant but the observations 
we wish to understand are made in this coordinate system. In Sections 4 and 5 we use the observed 
variation of power with b in RTN coordinates to show that on average the turbulent fluctuations 
are axi-symmetric and elliptically polarized. The polarization ellipse of the ensemble average of 
turbulent fluctuations is aligned along unit vector axes we define which are themselves aligned with 
respect to B. 

Finally we convert the observed power spectral tensor from RTN into magnetic-field-aligned 
coordinates and compare to the derived reduced power spectral tensor in this coordinate system. 
Combining the results in both coordinate systems allows us to demonstrate that the locally averaged 
magnetic field is an accurate representation of the axis of symmetry of the turbulence and therefore 
to plot the true reduced power anisotropy of the magnetic field. We show that the poloidal scalar 
function, which includes all pseudo-Alfvenic fluctuations, can be separated from the toroidal func- 
tion, which includes all shear Alfvenic fluctuations, in field-aligned coordinates and that they have 
different magnitudes and power anisotropy. We also show that the magnetic helicity is predomi- 
nantly in fluctuations with wave vectors near to the plane perpendicular to B. These properties 
represent important tests that any turbulence theory must satisfy and the results presented here are 
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Distance (AU) Latitude (°) \B\ (nT) \V\ (km/s) p (cm" 3 ) V A (km/s) ft 
2 - 2.28 79 - 74.8 1.5 ±0.3 780 ± 20 0.5 ±0.1 45 ± 8 1.6 ±0.6 

Table 1: Spacecraft location and average solar wind conditions for the 50 days of Ulysses data used 
in the analysis. 

important for all kinds of astrophysical turbulence: the solar wind, solar dynamo, and interstellar, 
galactic and intergalactic magnetic fields. 

2 Measuring the Reduced Power Tensor 

We use one-second resolution magnetic field data from the Ulysses spacecraft from days 200 to 249 
(inclusive) of 1995 when the spacecraft was in a continuous polar fast stream characteristic of high 
latitudes at solar minimum. The location of the spacecraft and average solar wind conditions for 
this period are summarised in Table [TJ 

We measure all components of the power spectral tensor of magnetic fluctu ations using a complex 



Morlet wavelet decomposition of t he time series reported in RTN coordinates (jHorbury et all 120081 : 



Podestal . 120091 : IWicks et all l20ld ) . Wavelet coefficients Wi are calculated using the inver s e Fourier 



trans form of the Fourier representation of the Morlet wavelet (jTorrence Sz Compol . Il998l : iPodestal . 

lj). 



2009) with the Fourier transform of the magnetic field B 

(9_ „ \ 1/2 poo 
=j±j J ^(wji-^Me-^'^e^dw (1) 

with St the time cadence of the data (1 second), s is the wavelet (time) scale which is varied to select 

different frequencies /, related by s = " 0+ ^~ t ~ U) -, ujq = 6, H(oj) is the Heaviside step function, and 
i,j run over R,T,N. The measured power spectral tensor as a function of frequency is then: 

P ij (f,t)=w i (f,t)w*(f,t) (2) 

Note that the wavelet amplitudes Wi contain phase information and are complex, making Pjj a 
Hermitian tensor. Each power measurement Pij(f,t), can be associated with the direction of the 
mean magnetic field B(t) calculated using the same averaging envelope at time t. The field direction 
is defined by the angles 9 B and ^ as shown in Figure |4j 

B = \B\cos9 B R + \B\s\n9 B cos(f) B t + \B\sm6 B sm4> B N (3) 

and we define: 

*=!§[ w 

Although RTN coordinates are more awkward theoretically than field-aligned coordinates for 
understanding magnetic turbulence, we use RTN to simplify data handling. The power contributions 
are accumulated and averaged in 404 separate direction bins: 18 equally wide 10 degree bins in 9 B , 
and variable width bins in <j) B to keep the solid angle area of each bin approximately constant. We 
keep bins equally spaced in 9 B since we are interested in the behaviour in this direction for physical 
reasons, thus where 9 B is near 0° or 180° there are fewer bins in <p B . The mean and standard error 
of each of the mean power contributions Pij = WiW* in each bin are then associated with the 9 B 
and 4> B at the center of the bin. Thus Pij(f, t) is converted into Pij(f, b) by this averaging process. 
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Figure 1: Real component of the power spectral tensor from Ulysses magnetic field data at / = 0.098 
Hz. Black areas represent bins that have fewer than 10 points in them. Red represents positive and 
blue negative contributions to the power, with white being zero. The colour scale has been scaled 
to the standard deviation of all points contributing to the power in each map individually, the value 
of which is shown above each panel. 



Using this method a map can be made of the power distributed over 9 b and 4>B for each of the 9 
tensor elements of Equation [2] at each wavelet scale. We have measured the anisotropy of the power 
spectral tensor at a range of frequencies (0.25 x 10~ 2 < / < 0.25 Hz), which allows us to verify the 
results presented here as typical over the inertial range of turbulence. In this paper we concentrate 
on a single scale since we are interested in power anisotropy; we will return to the scaling with / of 
the power spectral tensor in a future publication. Examples of the distribution of power in real and 
imaginary parts for each tensor element are shown in Figures Q] and [21 at a frequency of / = 0.098 
Hz, which is at the high frequ ency end of t h e anis otropic inertial range, maximizing the observable 
power anisotropy, as shown in IWicks et alJ l|201Gh . 

By definition the diagonal terms of the power tensor are real, since they are the wavelet coefficient 
multiplied by its complex conjugate. The off-diagonal terms are complex and have both real and 
imaginary parts. The data are presented in Figures Q] and [2] as a 2D map of the surface of a sphere, 
the horizontal direction in each of the nine plots is the <j)B direction and has variable bin- width, and 
the vertical direction is the 0b direction and has fixed bin width of 10 degrees. The color scale runs 
from dark blue for the largest negative values of power, through white at P = and then to dark 
red for the largest positive values of power. The colour scale is shown at the side of the plots and is 
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Figure 2: Imaginary component of the power spectral tensor from Ulysses magnetic field data at 
/ = 0.098 Hz. Black areas represent bins that have fewer than 10 points in them. Red represents 
positive and blue negative contributions to the power, with white being zero. The colour scale 
has been scaled to the standard deviation of all points contributing to the power in each map 
individually, the value of which is shown above each panel. Note that the diagonal components are 
zero by construction. 

calculated individually for each map in terms of the standard deviation a of all data contributing 
to that map (this includes the systematic and sinusoidal variations and so is considerably larger 
than the standard deviation in any individual bin, as shown in Figure [3|). The colour can be scaled 
onto the absolute value of the power by the value of a in nT 2 /Hz shown above each map. The off- 
diagonal maps appear noisier than those on the diagonal because their generally smaller magnitude 
makes the errors proportionally larger and they have both positive and negative regions with zero 
in between making any uncertainty in these regions appear more clearly in the color map. There 
are also fewer points per bin on average for 9 b > 90° making the error and therefore the scatter 
proportionally larger in this region. 

We will return to the absolute values of the data later, but for now we note that the maps have 
different, but clear, harmonic variations with 4>b- Within experimental error it appears that the 
RT and RN components are first harmonics of 4>b, their amplitudes are equal and they are 90° out 
of phase. The TT, NN and the real part of the TN components are second harmonics of <pB with 
similar amplitudes and multiples of 45° out of phase with each other. In Figure [2] the imaginary part 
of the TN element varies in only the Ob direction. Another striking feature of Figure [2] is that the 
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standard deviation of each map, used to calibrate the colour scale, is almost equal across all maps. 
The tensor is Hermitian by construction which provides the mirror symmetry about the diagonal. 
In fact, the dramatic dependence on 4>b is an artifact of using RTN coordinates, although it can be 
modified by the presence of non-axi-symmetric turbulence. In Section 3 we show how this arises, 
how it can be removed and how it can be used to extract extra information about the structure of 
the turbulence. 

We can quantify the (j>B dependence of the tensor components by fitting functions to the observed 
variation of power with cpg. We fit sinusoidal functions of cps to each tensor element containing 
a contribution from the R direction and sinusoidal functions of 2(f) b to the others, at each / and 
9b, with a non- linear least squares fitting method to determine the fitting parameters at each / 
and 6b as specified in Equations [5] - [TUJ The resemblance of these sinusoidal functions to the data 
provides motivation for employing them in fits. Their suitability is deeper than this, however, as 
we prove in Sections 3 and 4. The fits all consist of a real constant average A independent of 4>b 
and a real sinusoidal amplitude B as well as an imaginary constant average C and an imaginary 
sinusoidal amplitude D; each sinusoidal function also has a phase offset, E for the real part and 
F for the imaginary. This process is repeated for all values of 9b at which there are five or more 
points to fit to. As 9b coverage is limited by the reduction in solid angle close to 9b = 0, this means 
that the angle range covered is 15° < 9b < 175°, at each frequency, so there are never more fitting 
parameters (3) than data points (a minimum of 5), although the fits with 9b closest to zero are the 
least accurate. Thus each of the six independent elements of the tensor at each / and most 9b can 
be described by these scalar parameters with separate averages, amplitudes and phase shifts for the 
real and imaginary parts. 



Prr{(Pb) =Arr + B RR sin(</>£ + E RR ) (5) 

Ptt(4>b) =A T t + B T t cos(2^£ + E T t) (6) 

Pnn{4>b) =A nn + B NN cos(20 B + E NN ) (7) 

Prt{4>b) =A rt + B RT cosOb + E RT ) + i(C RT + D RT sin(</>£ + Fr T )) (8) 

Prn{4>b) =A rn + B RN s'm((pB + E RN ) + i(C RN + D RN cos((p B + F RN )) (9) 

Ptn(4»b) =A T n + B T n sin(20 B + E TN ) + i(C T N + D TN sin(2^ B + F TN )) (10) 



Tensor element A B CD 

RR 13.0 ±0.5 -0.1 ±0.8 

TT 15.8 ±0.8 -11 ±1 

NN 16.5 ±0.7 11.6 ±0.9 

RT 0.1 ±0.4 -2.4 ±0.5 0.03 ± 0.2 -1.0 ±0.3 

RN 0.3 ±0.4 -2.7 ±0.6 -0.05 ± 0.2 1.1 ±0.3 

TN 0.09 ±0.4 -11.0 ±0.6 2.2 ± 0.2 -0.1 ±0.3 

Table 2: Values for power in units of 10 -3 nT 2 Hz~ 1 of the fitted parameters in Equations l5l - [TUl at 
9b = 65° and / = 0.098 Hz, the same data as in Figure El The fitted angular phase is not shown 
since all are within errors of 0. Note that within errors there are only 4 independent real values and 
two independent imaginary values. 

Figure [3] shows a typical example of how the real and imaginary parts of each element of the 
power spectral tensor at a certain 9b, vary with 4>b- We can see that Re\PRr\ and Rc[Prm] 
elements behave to a very close approximation like cos<pB and sin^^ respectively implying that 
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Figure 3: The <ps dependence of power at 9b = 65° and / = 0.098 Hz. The blue and red points 
are the real and imaginary data respectively and the green and black lines are their sinusoidal fits 
as described in Equations [5] - [10] with the values quoted in Table [2j 



E ~ 0. Similarly the real part of Ptn, and Pnn and Ptt behave to a very close approximation like 
sin(20fi) and cos(20s) respectively, again implying that E ~ 0. We quantify this further by looking 
at the measured phase shifts E from the fitting of the real sinusoidal functions. The average shift 
of the phase in <ps over all scales and angles is only E = 0.6 ± 1.1° and always in the range ±10°. 
This is smaller than the angular resolution of the method we use (10°) and so within the accuracy 
of the method we cannot distinguish E from 0. 

Looking at Figure [3] again, we see that the imaginary part of the power varies to a very good 
approximation like ztsin^s in the Prt and Ptr elements and ± cos (f>B in Im[PfiN] and Im[PNR\- 
This implies that F ~ and again we quantify this by looking at the phase shifts F of the imaginary 
parts of Prt and Prn, which are too small to measure using our technique, being F = 0.2±1.0° and 
always in the range ±10°. This average ignores the Im[PrN] component since there is no sinusoidal 
variation and so these values of F are poorly constrained. 

Table[2]shows the fitted parameters for the data in Figure[3]at / = 0.098 Hz and 9b = 65°. There 
are only 4 measurably distinct non-zero real parameters (^4r_r, Att = Ann, P>tt = —Pnn = Ptn, 
and Prt = Prn) and two non-zero imaginary parameters (Drt = —Drn and Ctn)- This is true 
at all / and 9r in the range we studied. This means that there are at most six functions of 9b at each 
/ which together completely describe the properties of the power spectral tensor of the turbulence. 
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We show in the next section that this organization follows from geometry, the solenoidal field, and 
the conversion from power spectra in wave vector k to power spectra of the time series as seen at 
the spacecraft by deriving six (/^-independent functions corresponding to the six ^-independent 
values observed in the data. 



3 Systematic Effects of Geometry 



There are several potential problems in comparing observations with theory in solar wind turbulence 
studies. One is that all in- situ spacecraft observations of the solar wind are of a 'reduced' spectrum, 
but theory usually addresses the spectrum in k space. Power spectra calculated from time series 
of single point observations made in a fast flowing medium carrying a relatively slowly evolving 
turbulence are, by Taylor's hypothesis ( Tavlorl . ll938l ). an integral of the Pjj (k) in wave vector space 
over the plane perpendicular to the flow defined by k ■ V = 2irf ( Fredricks fc Coroniti . 19761 ): 



Pij(f,V) 



Pij{k)5{2nf -k-V)d 3 k. 



(11) 



When Pij (k) is anisotropic, Pij (/, V) will depend on the direction of V relative to any symme- 
try in Pij(k). If symmetry in Pij(k) is organized by the direction of the local mean magnetic 
field, b, Pij(f,V) can be better written as Pij(f,b). In fact the total (trace) power P(f,b) = 
y^! Pji(.f, b)„-_ RT A T is known to be anisotropic in both power and spectral i ndex as a function o f 9 b 
(jHorburv et all . I2OO8I : IPodestal . 120091 : ILuo fc Wul . I2OI0I : IWicks et all I2OI0I : iForman et all . l201lh . In 



order to derive the expected geometrical effect of the reduction on the spectrum we must account for 
three vectors (k,B, V) and their corresponding coordinate systems aligned with b and V. Figure 
H] shows these vectors, V the solar wind velocity, b the unit vector of the magnetic field, both of 
which are measured in RTN coordinates in this analysis, and k the unit wave vector of a turbulent 
fluctuation. 




Figure 4: The geometry used in this paper, V is the solar wind velocity, considered to be in the 
radial (R) direction, b is the unit vector of the magnetic field, k is the unit wave vector of a 
fluctuation. R, T and N are heliocentric coordinates, t and p are toroidal and poloidal directions, 
9b and 4>b are angular coordinates of 6, and 9^ and (j)^ are the angles between k and the (V,b) 
plane. 

The magnetic field is solenoidal, V • B = 0, ensuring that all fluctuations are confined to the 
plane perpendicular to k, which is tangent to the surface of the sphere at position k in Figure [H 
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To describe the fluctuations of SB in this plane we define the toroidal direction i, perpendicular to 
both k and b, and the poloidal direction p, perpendicular to k and t thus: 



b x k 



b x k 
k xi 



(12) 



(13) 



Considering the sphere in polar coordinates with b as the polar axis and k as the radius vector: p is 
in the direction of decreasing 0^, and t is in the direction of increasing 0^; (fc, t,p) is a right-handed 
coordinate system. 

The power spectral tensor of toroidal fluctuations alone is a scalar function Tor(k) times the 
dyadic [t : t], and for poloidal fluctuations alone is [p : p] Pol(k). If both polarizations exist, any 
correlation between them will result in the additional power spectral elements It : p-\-p ■ t] C(k) 
which is real and i\t:p — p:t\ kH(k) which is imaginary and anti-symmetric. The four scalar 
functions are all real and the spectral tensor structure correspon ds in detail to t he co mplete 
description of transverse light waves with the Stokes parameters ( Chandrasekhar . 196dl ) (/ = 
Tor + Pol, Q = Tor - Pol ' U = 2C \V = 2kH), and to the description of solenoidal MHD fluc- 
tuations by Oughton et al.l ( 19971 ) where our Tor(k), Pol(k), C(k) and H(k) correspond to E, 
E - (6 x k) 2 F, k{b x kfC and H in their paper. The toroidal fluctuations are perpendicular 
to b, and so they are sometimes called 'Alfvenic', since this is the polarization of small-amplitude 
shear Alfven waves. Similarly, the Pol(k) fluctuations are sometimes called 'pseudo-A lfvenic' be- 



cause their polarization is the same as that of small amplitude pseudo- Alfven waves (c.f. ICho et al 
20021 ). Following these theoretical structures we define the anisotropic power tensor as a function 
of k: 

P{k) = Tor(k) [t:i]+ Pol(k) [p : p] + C(Jfe) [i : p + p : i] + ikH(k) [i : p - p : i) . (14) 

This formalism is completely general and Equation [T5] describes any turbulent field satisfying the 
solenoidal condition regardless of any symmetry. 

In order to use Equations (fTTI - [T4|) with the measured tensor in RTN coordinates we must 
express t and p in RTN coordinates, but keep k in field-aligned coordinates. For this we define a 
new coordinate system aligned with b and containing the radial flow direction of the solar wind: 

e z = b (15) 



e z x V 
\e z x V\ 



x R 



sin I 



k 



(16) 

(17) 
(18) 



Defining e y this way means V is in the x-z plane and thus V k = \V\ (s'm0Bk x + cos#b&z)- We will 
exploit the fact that k y is therefore not in the delta function of Equation [11] to help us understand 
the symmetries of the scalar functions later. It should also be noted that <f>s is not in the delta 
function either, thus <pB can potentially be moved outside the integral in Equation II li 

We find t and p in RTN coordinates for any k using Equations [15] - [18] in Equations 1121 and [TBI 
with k i = 



k^. ~\~ ky and k 



U2 _i_ Ul _i_ VI. 
x ' y ' z ' 



( 



-%-sm0 B R 



\ 



sin <pB + cos 9b cos 

' B + 1?-COs6b SU10B 



COS i 



T 

N 



(19) 
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p 



sin 9 B cos 4> B - jjTj- (k y sin 4> B - k x cos 9 B cos <f> B j) T 
, sin sin 4>B + ^ (ky cos B + fca cos 6> B sin 4> B )j N 



(20) 



Note that <p B does not appear in the R component of either t or p, and only as sin^ or cos<ft B 
in the other two components. Since <p B is not involved in the integration in Equation QTJ this 
dependence appears directly in the maps of P%j(f, b) as 0, 1 st , or 2 nd order harmonics of (j) B , and is 
easily seen in the maps of Pij(f, 6) in Figures [JJ and [2] The amplitude of each of these harmonics 
in <j) B is an integral involving 9 B over the power distribution P(k) and is a function of 9 B . We will 
now determine what these amplitude functions are by gathering terms with no dependence, first 
harmonic and second harmonic dependence on (f> B . 



4 Harmonics of 



'B 



There are only six independent linear combinations of the elements of P RTN (f,b) which have no, 
first or second harmonic 4> B dependence. This is exactly the same number as the independent power 
amplitudes of the fitted sinusoidal functions of <f> B observed in Table [2j These can be expressed in 
terms of the four scalar functions by putting Equation [14] through the reduction integral (Equation 
TT]) with the t and p vectors in RTN coordinates, as in Equations [19] and [20] The results are 
projections of the four scalar functions from fc-space onto the RTN coordinate system which we 
now derive. 

The Trace and Prr are independent of <p B and Equation 1111 gives the two projection integrals: 
PrrUJb) = 

f t R Tor(k) + p 2 R Pol(k) + 2t R p R C{k))8{2TTf — k- V) d 3 k, (21) 



Trace(/,fl B )= ]T Pu(f,b) 



i=RTN 

Tor(k) + Pol(kj\ 5(2vr/ -k-V) d 3 fc, (22) 
We then collect the real components which are first harmonics of <p B , if we combine them we find: 
h e^ B = Re (p RT (f, 6)) + iRe (P RN (f, 6)) = 

t R Z t Tor(k) + p R Z p Pol(k) + (t R Z p + PR Z t ) C(k)) 6(2nf -k-V) d 3 k (23) 




where we have expressed the <f> B dependence of t and p compactly in the following way: 



Z t (9 B ,k)e t4,B = t T + it N = e^ B -^{k y cos 9 B - ik x ) (24) 

k± 

Zp{9 B ,k)e l * B = p T + i PN = e **B (ft s [ n 9 B + k x k z cos 9 B + ik y k z ) (25) 



Note that Zt, Z p are complex functions of 9 B and k so I\ is complex and the absolute phase of P R t 
and P R n are important clues to finding Tor(k), Pol(k) and C. 
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Similarly collecting the real components which are second harmonic in ^ we find: 

j 2 e 2i<t>B = p Tr(/) b) _ p NN ( fi h) + % (p TN {f, b) + p NT (f, i>\ 

Z?Tor(k) + Z*Pol(k) + 2Z p Z t C(k))s(2irf - k ■ V)d 3 k. 




(26) 



this integral is also complex and so the absolute phase is also important. 

From the imaginary part of P RTN (/, 6) we find there are only two equations with simple de- 
pendence on <f>B, both projections of the scalar function H(k). The tensor form multiplying H(k) 
in Equation 1141 is: 

1*iPj Pi^j — €ijm\P ^ P)m — Cij'm^m (^7) 

so that 



Im 



Pii(f,b) 



k m H(k)5(2irf - k ■ V)d 3 k. 



(28) 



Combining the two off-diagonal terms with first harmonic dependence on (fig: 



h e 



l<pB 



Im[P RT (f, b)] + Um[P RN (f,b)} 



(t R Z p - PR Z t ) kH(k)5(2irf - k ■ V)d 3 k 



(29) 



this integral is potentially useful as it can give us information about the symmetries of H(k) through 
the different projections of t and p it contains. 

Finally, we find that one imaginary term has no dependence on </>#: 



Im 



PtnV, b) 



(k ■ R)H(k)5(2irf - k ■ V)d 3 k, 



(30) 



a nd by noticing tha t k ■ V = \ V\(k ■ R) we easily recover the well known and frequently used result 
of iMatthaeus et al.1 (Il982l ) that Ii7i[Ptn] is the reduced magnetic helicity: 



Im 



PtnU, b) 



2vr/ 
V 



H(k)5{2irf - k ■ V)d 3 k. 



(31) 



We have not yet made any assumptions about the symmetry or behaviour of the four scalar 
functions and so these six relations are generally true. We now look at the Ulysses data in more 
detail to see what restrictions the observed tensor elements place on the four scalar functions. 



5 Observational restrictions on the scalar fields 

We can now return to the observations in Figures Q313] and Table [2] and compare them to the geo- 
metrical effects derived above. We can rearrange Equations [23] and [26] to illustrate the dependence 
of the individual tensor elements on 4>b and the integrals I\ and I2 which are functions of 9 R and 
/. We drop the dependence on / and 9b temporarily for simplicity: 



Re[P R x\ = Re[Ii] cos{4>b) — Im[I\] sm{4> 

Rc[Prn] = Re[h] sm((pB) + Im[Ii] cos((p 

Re[P NT ] = ^(i?e[/ 2 ]sin(2^ B ) +/m[/ 2 ]cos(20 B )) 

Re[P TT - P NN ] = Re[I 2 ] cos(20 B ) - /m[/ 2 ] sin(20 B ) 



>b) (32) 
>b) (33) 

(34) 

(35) 
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A similar rearrangement of the imaginary part using Equation 1291 yields: 



Im[P RT ] = Re[h] cos(4> B ) - Im[h] sin(^ s ) 
Im[P RN ] = Re[I 3 ] sinOs) + Im[I 3 ] cos(<^b) 



(36) 
(37) 



A surprising property of the data described earlier is that the angular phase offsets E ~ 
and F ~ for all of the measured quantities. By comparing Equations [5] - [10] with E set to 
with Equations [32] and [33] we can see that in the solar wind this implies that ~ and in 

addition using Equations [33] and [35] Im[I<2\ ~ 0. F ~ and Equations [36] an d [37] similarly imply 
that Re[I 3 ] ~ 0. The implications of this are discussed later in this section. 
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Figure 5: The amplitude of the fitting parameters A, B, C, and D for the six independent power 
spectral tensor elements as a function of 9b, at / = 0.098 Hz. 



Equations [38] - 03] show the resulting 6b and (f>B dependencies of each element of the power 
tensor measured in RTN coordinates on the six ^-independent integrals defined in the previous 
section. These functions fit the observations very well as can be seen by comparing the equations to 
Figures [H [2] and [3] Just as in Table [2] there are only four real and two imaginary amplitudes. These 
terms are a result of the geometry shown in Figure H] and the reduced nature of the measurements 
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as described by Equation [TT] as well as the turbulent power spectrum. 



Prr =Prr(0b) 
1 



Ptt =2 [T±bxx(0b) ~ Prr(0b) + h(0 B ) cos(2<^ B ) 
Pnn =^(Trace(^ B ) - P RR (6 B ) - I 2 (8 B ) cos(2<£ B ) 

^Prt =I\{0b) cos(0 b ) - ih(0B) sin(0 B ) 
fRJV =I\{0b) sin(0 B ) + ih{0 B ) cos(0 B ) 

i 3 ™ =V/ 2 (^)sin(20 B ) + iff m (0 B ; 



(38) 
(39) 

(40) 

(41) 
(42) 

(43) 

We can now go back to the results and look at the amplitude of the fitting parameters A, B, C, 
and D (Equations [5]- [10] and Table [2]), as a function of 9 B , shown in Figure [5] The error bars are the 
standard error from the linear least squares fitting. Figure [5] shows a remarkable amount of variety 
in the power anisotropy of the different tensor elements with 9 B , including that all of the diagonal 
terms have less power in the field para ll el dir ection (0r — > 0°) than in th e perpendicular direction, 
recovering the results of Bieber et al. ( 19961 ) and Horburv et al. ( 20081 ) that the total power is 
anisotropic. We can see again that there are only 4 measurably distinct non-zero real parameters 
A RR , Att = Ann, Ptt = —Pnn = Ptn, and B R t = B R n and two non-zero imaginary parameters 
Drt = —Drn and Ctn and that their dependence on 8 B is similar. Rather than interpret this 
in terms of P RTN (f,6 B ) the direct link with the projections of the four scalar functions is more 
obvious if we consider the six ^-independent functions in Equations EH [22] [23] [26] [29] and [3D 

We can extract the six ^-independent functions from the fits to the data in Section 2, and can 
even measure /1-3 in two different ways to cross check the results: 



Prr{0 b 
Trace(8 B 
H m (9 B 
h(0 B 
h(0 B 

h(0B 



--A RR {9 B ) 

--A RR {6 B ) + A TT {6 B ) + A NN (9 B ) 
-2Ctn(0b) 

-B R t{0 b ) = B R n{6 b ) 

-Btt(0b) — Pnn(0b) = 2Ptn(6b) 
- - D R t(8 b ) = D R n(0 b ) 



(44) 
(45) 
(46) 
(47) 
(48) 
(49) 



Figure [6] shows these six ^^-independent functions directly measured from the fitting of the power 
spectral tensor. The two different ways of measuring /i_3 all agree with each other remarkably 
well, which is not re quired in general, but shows that the turbulence is solenoidal and that the 
Ouehton et~"all (1 1 9971 ) theory applies and our subsequent derivations are correct. It is also interesting 
to note that although I\ and ^3 appear sinusoidal in 28 B upon closer inspection the peaks of the 
power are shifted from 45°. 

By considering the properties of Trace, P RR , I\ and I2 we can deduce some important properties 
of the underlying Tor(k), Pol(k) and C(k) functions. As discussed above we observe that within 
errors I\ has no imaginary part. We also know that Tor(k) and Pol{k) are not both identically zero 
since both the Trace and P RR are not zero. If we write out the complex terms in Equation [23] for 
I\ we find that the imaginary part contains terms in Tor(k), Pol(k) and C(k) all with pre-factors 
that are odd functions of k y . One possible way for Im[Ii] = is therefore if Tor(k) and Pol(k) are 
certainly even functions of k y , that is mirror-symmetric about the (V,B) plane, sin ce they then 
integrate to zero. Tor(k) and Pol{k) are shown to be even in lOughton et al. I (|l997h and so our 
results are in accord with theirs. Axi-symmetry about b is a stronger conclusion not proven, but 
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consistent with mirror-symmetr y in k y . Similar l y C(k ) must be even in k y , or alternatively it can 



be zero. However, as shown in lOughton et al.l (|1997l ). C(k) is necessarily odd, so that we must 
conclude C(k) = 0. 

Applying the same analysis to the imaginary part of I2 we see that the terms multiplying Tor(k) 
and Pol(k) are also odd functions of k y and so are integrated to zero by Equation II 1} however the 
terms that multiply C(k) are even in k y and so mirror-symmetry cannot be used to explain the lack 
of contribution to the power. Thus since Im[I\{f, 9b)] = Im[l2{f , 9 b)} = we again conclude that 
C(k) = 0. 

Continuing to ^3 we find that the real part, which involves only the magnetic helicity H(k), has 
pre-factors that are odd functions of k y , and the imaginary part has even pre- factors in k y . We 
observe that itefis] = since F ~ 0, so again following the same line of reasoning we find that 
H(k) is also even in k y and could be axi-symmetric about b. 

This is as far as we can conveniently proceed with the analysis in RTN coordinates. We have 
used the geometrically induced (fis dependence to draw conclusions about the symmetries of Tor(k), 
Pol(k) and H(k) and to show that C(k) = 0. The RTN coordinate system, however, also imposes 
strong geometrical 9b dependencies, as can be seen in Figures [5] and [H The large imposed 9b depen- 
dence of Prr, and /1-3 (Equations \21 \ [23 l 1261 and I29 p make it hard to extract further information 
about Tor(k), Pol(k) and H(k) from them. Therefore, at this point we convert the observed power 
spectral tensor in RTN, to the ^-aligned XYZ coordinates. 
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Figure 7: The real component of the power spectral tensor from Ulysses magnetic field data at 
/ = 0.098 Hz in field-aligned-coordinates. Black areas represent bins that have fewer than 10 points 
in them. Red represents positive and blue negative contributions to the power, with white being 
zero. The color scale has been scaled to the standard deviation of power in each map individually, 
the value of which is shown above each panel. Compared to Figure Q] in RTN coordinates much of 
the variation has been removed. 

In the field-aligned XYZ coordinates defined in Equations 1151 - [T71 the t and p vectors are: 



i - -h-e. +—e 

~ k ± x+ k ± y 

k x k z kyk z k_\_ 

p ~ ~Jk^ ex ~Jk^ ey + T e 



(50) 
(51) 



importantly they do not involve 9b or 4>B 
to XYZ coordinates is 



at all. The transformation of a wavelet coefficient in RTN 



ijWj 



(52) 



and 



where j = R, T, N, and i = x, y, z 

sin#B — cos 9b cos 4>b — cos 9b sin 4>b 

M= | sin 0b — cos 4>b 

cos 9b sin 9b cos (f>B sin 9b sin (pB 



(53) 
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We use this matrix to transform the observed reduced power spectral tensor from the RTN coordi- 
nates of the data to the XYZ coordinates aligned with b below, using: 



P xyz (f, B ) = MP RTN (f, b)M T . (54) 

In field-aligned coordinates the transformed power spectral tensor should have no 4>b dependence 
and simpler 9b dependence than shown in Figures [I]|3] since the t and p vectors in XYZ are no 
longer dependent on these angles. 

The expressions for P xyz {k) using the definitions of t and p in Equations 1501 and 1511 are: 

P xx {k) = ^Tor(k) + %Pol{k) + 2 k -^C(k) (55) 

P yy (k) = ^Tor(k) + %Pol{k) - 2 k ~^C(k) (56) 

P zz {k) = ^Pol(k) (57) 



k x k y _, , k x k v k z t ^ (k y k x )k z 
_^Tor(fc) + -^-- 



P*vW = '^Tor{k) + ^^Pol(k) + 11 C ik) + ik z H(k) (58) 



P X z{k) = - -j-fPol(k) - -fC{k) - ik y H(k) (59) 
p y^ k ) = ~ ^T Pol &) + jrC(k) + ik x H(k) (60) 

Although there is explicit dependence on components of k in every term there are no 9b or <pB 
dependent geometrical pre-factors. C(k) and terms odd in k y are shown for completeness, although 
we know from the RTN analysis that they do not contribute to P x v z (f : 9b) in the solar wind. The 
reduction to P(f,9s) using Equation [TT] introduces dependence on 9b if there is any anisotropy in 
the scalar functions. 

Figure [7] shows the same results as Figure [H the real part of the power spectral tensor, converted 
to XYZ coordinates. There is no (j>B dependence in any of the elements and no sinusoidal-like 
dependence on 9b- This result, combined with the strong 4>b and 9b dependence of P RTN implies 
that we have correctly identified the direction of b using the local wavelet averaging method. If we 
had identified b incorrectly with a systematic error then there would be a sinusoidal dependence 
on 9b in P x u z and if there was a random error smoothing out the variations we could not have 
measured the precisely predicted <j>B dependence of P RTN . 

Figure [8] shows the 9b dependence of each independent element in both the real and imaginary 
parts; the error bars are calculated as the error on the mean of all data contributing to each bin in 9b- 
This figure allows us to make further deductions about the scalar functions using a similar analysis 
procedure to that of the RTN tensor previously. First, the real parts of the off-diagonal elements 
are all observed to be within errors of zero. Looking at Equations [58] - EQ] we see that Tor{k) only 
appears combined with an odd function of k y and so we have rediscovered its mirror-symmetry, 
Pol(k) is also combined with an odd function of k y in two of the elements but P xz (k) oc k x k z Pol(k) 
so Pol(k) must be predominantly 2D, that is mostly confined to \k z \ ~ 0. Finally C(k) is multiplied 
by k y and k y in P xz {k) and P xy {k) respectively and so must be zero, as previously discovered. 

The imaginary parts Im[P xy ] and Im[P xz ] are within errors of zero. Im[P xz ] is combined with an 
odd function of k y so again the data imply that H(k) is mirror-symmetric in k y . Im[P xy ] = implies 
a further symmetry of H(k) in k z , similar to Pol(k) this implies that in the inertial range H(k) 
is mostly associated with fcs in the plane perpendicular to B. At the highest frequencies studied 
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Figure 8: Power anisotropy of each power spectral tensor element as a function of 9b in magnetic- 
field- aligned coordinates at / = 0.098 Hz. The error bars are calculated as the error on the mean 
in each 9 b bin which have no dependence on 



here Im[P xy ] becomes anisotropic (a very small signature may be visible in Figure [8] although it is 
indistinguishable from when errors are considered and is zero in the inertial ra nge in general) 



however this is p r obabl y associated with plasma instabilities near the ion gyroscale (|He et al.l . 12011 
Podesta Sz Gary . 201 ll ). Only Im[P yz ] oc k x H(k) has finite value in general, although it is much 



smaller than the Trace and has much weaker 9b dependence than IttiIPtn]- 

Moving on to the diagonal elements, we see that as in RTN coordinates the diagonal elements 
are anisotropic with more power at 9b ~ 90° than at 9b ~ 0°. They are ordered in power with 
Pyy > P xx > P zz seeming to show a 3D anisotropy with the most power perpendicular to the (V, B) 
plane, then intermediate power perpendicul ar to B but in the (V , B) plane and finally the least 
power parallel to B, exactly as first shown by lBelcher &: Davisl ()197ll ) and in agree ment with the well 
established results that the solar wind is an i sotropic with < P 7Z < P TiT +P ,„, (e.g. iMatthaeus et al 



199d ; basso et all 120051 : IMatthaeus et all 120051 : lOsman fc HorburyLl2007l ). Looking at Equat 



10ns 



]we see that this ordering implies that < 



%Pol{k) < Tor{k) + %Pol{k). 



One way this can 



be achieved is if Pol(k) < Tor(k), although this is not required. 

In the assumption that the turbule nce is 2D the reduc e d power P vv (f,0Fi) > P xx (J,9b) is due 
to the reduction integral (Equation [TT]) ( Bieber et al.l . Il996l : iTurner et al.l . 1201 ll ) . We can now make 
a stronger statement than in the previous paragraph, we have shown that Pol(k) is mostly due to 
approximately 2D wavevectors and so does not contribute strongly to either of these terms since they 
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both contain k^Pol(k), thus the observed power is from the purely Alfvenic Tor(k) fluctuations. 
Pxxifi^B) and P yy (f,0B) are reduced in Equation [TT1 as oc kyTor(k) and oc k^.Tor(k) respectively 
and it is inter esting to note that P yy ( f, Ob)/ Pxx{f, Ob) is approximately 2 at all values of Ob at this 



frequency. As lTurner et al.l (|201ll ) showed the reduction integral applied to a power spectrum of 2D 
fluctuations results in a constant factor proportional to the spectral index of the turbulence, so this 
apparent anis otropy is a featur e asso ci ated with sampling a long a single cut thr ough the data and 



the results of lBelcher fc David (|l97lh : bieber et all (|1996h : iTurner et all (l201lh now have unified 



explanation. We have also separated the reduced form of the Alfvenic Tor{k) with two different 
projections from the reduced pseudo- Alfvenic Pol(k) in the observations. They have different power 
levels and anisotropy since P zz (f, Ob) is very different in magnitude and shape from P yy (f, Ob), which 
is approximately 2P xx (f,0 B ). 

6 Summary and Conclusions 

We have observed all nine elements of the reduced power spectral tensor of MHD-scale fluctuations 
in fast solar wind using wavelet transforms of magnetic field observations by the Ulysses spacecraft. 
Each element of the tensor is resolved using angle coordinates Ob and ^ at a single frequency 
/ = 0.098 Hz. The signal is anisotropic and depends on the direction of the local mean magnetic 
field. This anisotropy can be seen in Figures Q] and [2] and it is quantified in the 4>b direction by 
fitting sinusoidal functions in Figure [3] with the amplitudes given in Tabled These show that within 
errors we observe only six ^-independent power amplitudes in the solar wind. 

We explain the generation of thi s anisotropy analytica lly by applying a scalar field and tensor 
description of solenoidal turbulence ( Oughton et all 1997I ). We choose our scalar functions so that 



they represent the toroidal (Tor(k)) and poloidal (Pol(k)) fluctuations with respect to the local 
mean magnetic field direction b, and their in- and out-of-phase correlations (C(k) and H(k)). We 
then convert this representation into the spacecraft data coordinate system where the toroidal t 
and poloidal p directions are expressed in terms of heliocentric RTN coordinates (Equations 1191 and 
20|) . Applying the reduction integral to the four scalar fields in conjunction with the appropriate 
dyadics of t and p we derive the dependence of the reduced power tensor Pij(f, b) on the four scalar 
fields. 

While we do not know the analytical form of any of the four scalar fields, the geometrical 
dependence on in the RTN coordinate system is independent of the reduction integral and 
simple sinusoidal dependences on 4>b are found. To simplify the equations we gathered terms with 
no, first harmonic or second harmonic dependence on (j) B . We found six combinations of the reduced 
power spectral tensor that analytically have precise sinusoidal (j>B dependence, matching up with 
the six observed independent power amplitudes. These are the Trace, Prr and IthIPtn] elements 
of the reduced power spectral tensor which do not depend on 4>b, and the I\, ^and 1% amplitudes of 
the harmonic components in the other elements, defined in section 4. Observations of these from the 
solar wind do indeed have no (fts dependence confirming the derivations and their Ob dependence 
is shown in Figure [6J From these results we draw several conclusions: 

1. The 4>b dependence of the data (Figure[3]) follows the form of the power spectral tensor derived 
for solenoidal fluctuations transformed from field aligned coordinates in to RTN coordinates. 

2. Im[I\] = Im[l2\ = implies that the turbulent power in Tor(k) and Pol(k) is an even 
function of k y and so is mirror-symmetric about the (V, B) plane, and thus it is likely to be 
axi-symmetric about b. 
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3. The integrand of Im[Ii] oc k y C and the integrand of Im[l2\ oc k y C and both are observed to 
be zero when integrated; thus the scalar function C(k) = 0. 

The strong 9 b dependence observed in all panels of Figure [6] arises from the integrals over 



a combination of any actu a l anisotropy of the turbulence (IBieber et al.l. Il996l; iDasso et al.1. 120051 ; 



Goldreich fc Sridharl Il99,l Il997l : iMatthaeus et al 



199fil . Il998l : lOughton et all . Il998l . 1201 j ) and 

and [29] as their dependence on the scalar 



geometrical effects. This can be seen in Equations | 
fields, the coordinate transformed unit vectors t and p, and the complex functions Z-t and Z p 
(Equations 1241 [25]) . Since we cannot make simplifying assumptions such as symmetries in the 9 b 
direction, and since we do not know the analytical form of the scalar functions, further progress in 
the RTN coordinate system is difficult. 

We therefore transform the observed P(f, b) into magnetic-field-aligned coordinates XYZ. All 
the (pB variation in power disappears in agreement with the theoretical prediction (Equations 1551 
60|) . The observations have completely reproduced the theoretical prediction for the (ps dependence 
in RTN and the independence in XYZ coordinates, of the power, confirming that our measurement 
of b using the local mean magnetic field is an axis of symmetry for the ensemble average. This is a 
strong justification for using the local mean field when studying anisotropy in turbulence since our 
results indicate that this direction has a strong influence on the symmetry of the scalar functions. 

By comparing the analytically derived field-aligned power spectral tensor elements (Equations 
55ll60p with the data in Figures [7] and [8] we can draw further conclusions: 

4. The Alfvenic Tor(k) fluctuations can be separated from the pseudo-Alfvenic Pol{k) fluctua- 
tions since P zz (f,9s) is a function of Pol(k) alone. 

5. Pol(k) is measurable and has a different power anisotropy with respect to 9b than Tor(k). 
Since observationally P zz < P xx ,P yy it also seems likely that Pol{k) < Tor{k). 

6. Pol(k) is even in k y and therefore it is likely to be axi-symmetric about b since the observed 
reduced real parts of P X y{f,&B) and P yz (f,0B) are zero and P xy (k) oc P yz {k) oc k y Pol{k). 



7. Pol(k) is mostly due to fluctuations with \k z 
zero and P xz (k) oc k z Pol(k). 



because the observed reduced P xz (f, 9b) is 



8. C{k) = is confirmed by all real off-diagonal elements being zero since combined they have 
both odd (P xz and P yz ) and even [P xy ) pre- factors in k y and k x . 

9. H[k) is even in k y and therefore likely axi-symmetric about b, since the imaginary off-diagonal 
element P xz {f,9s) is zero and P xz {k) oc k y . 

10. H(k) comes from fluctuations that have wave vectors in the plane perpendicular to B with 
power confined around \k z \ ~ since the imaginary part of the off-diagonal element P xy is 
zero and P xy {k) oc k z 

A physical interpretation of these results is that turbulence in the solar wind is made up of mostly 
toroidal fluctuations that are anisotropic. The observed P zz and therefore Pol{k) is compatible with 
solenoidal fluctuations as in Equation 1571 however, a spectrum of \B\ fluctuations is observed in the 
fast solar wind, so we cannot rule out compressible plasma fluctuations as a source of this variation. 
If we consider the results in terms of a superposition of pola rized fluctua t ions then Tor(k) > Pol(k ) 
implies the fluctuations must be elliptical on average (e.g. IChen et all l201ll : iMallet et all Ell). 
H{k) ^ 0, implied by the finite values of PtnU ,&b) and P yz (j ,9b), means that the Tor(k) and 
Pol(k) fluctuations are partially correlated and there is a polarization ellipse. C(k) = implies that 
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the ensemble averaged polarization ellipse is oriented along i or p ( Chandrasekhar , 196dl ) but since 
the solar wind is not entirely coherent waves this must be a result of a superposition of polarization 
ellipse orientations that average to zero. Th us the ensemble average turbulence is similar to partially 
polarized, p artly natural ( incoh erent) light ( Chandrasekhar! . Il960l ). 

Recently iTurner et al.1 (|201ll ) showed that the difference in power P yy {f) > P X x{f) can arise from 
the reduction of an axi-symmetric 2D turbulence. Here we have shown why they find agreement 
between a superposition of 2D Alfven waves, numerical MHD simulations and the solar wind: the 
Alfvenic Tor{k) dominates the pseudo-Alfvenic Pol{k) contribution to both P xx and P yy when they 
are reduced and so observations of these terms appear A lfvenic, even if pseud o-Alfvenic fluctuations 
exist. The results we have shown here set the work of Turner et al. ( 201ll ) in the wider physical 
context of the full turbulent power spectral tensor. 

The observation process demonstrated in this paper can be repeated at many different scales 
rather than just one so the scaling of the <j>B invariant functions and the field aligned power spectral 
tensor can be measured. This may help test different theories for anisotropic turbulence if theoretical 
predictions for the scaling of the scalar functions are made and we intend to present such an analysis 
in the near future. Finally, this work also demonstrates that care should be taken when using off- 
diagonal ter ms from the pow er sp ectral tensor to observe physical phenomena. For example, in 
work such as iHe et al.1 (j201lh and IPodesta fc Garvl (j201lh . the m agnetic helicity is measured as 



Im \Ptn] /Trace, however the trace has it s own power anisotropy ( Horbury et al. . 20081 : Podesta, 



20091 : iLuo fe~Wul . l20ld : IWicks et all l20ld ) which depends mostly on Tor(k) which we have shown 
is different from the anisotropy of H(k) alone. Thus dividing by the trace introduces or removes 
apparent anisotropy from these results. Furthermore from the work presented here we can see that 
Pxyi Pxz and P yz contain different projections of H(k), from which we may learn more about the 
properties and symmetries of the helicity induced by solenoidal turbulence and instabilities. 
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